Mimicking Friedmann-Robinson- Walker universes with tunable cold Fermi atoms: 
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In a series of papers to demonstrate emergent 'gravity' some authors have explored phonon pro- 
duction in BECs with a tunable speed of sound, in particular to emulate FRW universes. The 
premiss is that, within such systems, the gapless mode (the phonon) looks Lorentzian for low mo- 
mentum at least. However, when it comes to phonon production in cold Fermi gases whose speed 
of sound is controlled by a Feshbach resonance, it is impossible to shake off the the underlying 
Galilean invariance because of the interplay between gapless and gapped modes. Such phonons as 
are produced do not follow the pattern anticipated for FRW metrics, at variance with the aims of 
the programme. 
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I. INTRODUCTION 

In this paper we analyse a Bose condensate as its speed 
of sound is manipulated experimentally. The reason for 
so doing is that the low-frequency long- wavelength modes 
of the phonon field 9 satisfy the geodesic equation 

g» v d v 6) = 0, (1) 



for some acoustic metric g determined by the local speed 
of sound. In consequence, the condensate looks to be 
an ideal candidate with which to study 'analogue grav- 
ity' or 'emergent spacetime', a programme begun more 
than thirty years ago, when Unruh observed [l[ that flu- 
ids moving faster than their local speed of sound gen- 
erated event horizons with many similarities to horizons 
in general relativity. [A recent review article on this 
programme has over 700 references.] 

The tactics for such an approach with condensates have 
been addressed by several authors [3141^. The speed of 
sound can be varied by changing the density of the con- 
densate (To} , but in this paper we examine that class 
of condensates in which the speed of sound is tuned 
through the application of an external magnetic field 
which changes the binding energy of a Feshbach reso- 
nance. This allows us to enter a strong coupling regime 
where the scattering length is larger than the interparti- 
cle separation. We are not interested here in event hori- 
zons jlOj but in the phonon production that arises as we 
change the metric, mimicking particle production in the 
early universe. Spatially homogeneous condensates with 
tunable sound speeds and, thereby, simple metrics g(t), 
look prime candidates with which to pursue analogies 
with FRW universes fllt-fli]. 

This paper is, in lar ge p art, a commentary on these 
latter three papers |12| - [l4| which estimate phonon pro- 
duction for such metric behaviour. As discussed in them, 
some caution is necessary. The equation (fT]) only de- 
scribes the linear long-wavclcngth part of the dispersion 



relation. In practice, non-linear dispersion sets in quickly 
at shorter wavelengths, since it is no more than the re- 
assertion of Galilean invariance of the condensate, tem- 
porarily hidden at long wavelengths. Nonetheless, it was 
anticipated in [l3|, [TJ] that there would still remain a nar- 
row regime in which the analogy with FRW universes is 
sufficiently good to mimic the creation of quantum modes 
in an inflationary regime. 

In order to make a comparison with (fT]), a robust 
analytic semi-classical approximation is necessary as a 
first step. Specifically, the analysis in fl2l - fl4( is predi- 
cated on the representation of the condensate by a Gross- 
Pitaevskii (GP) mean bosonic field in which the strength 
of the self-coupling, and hence the scattering length, is 
given an explicit time-dependence. In general, for cold 
bosonic atoms, increasing the scattering length increases 
the effect of three-body interactions [15j . As a result, 
even for dilute gases, it is difficult to find circumstances 
in which the GP mean field approximation is reliable (l5j 
in its unadulterated form. This is not a blanket prohi- 
bition. There are situations in which mean-field theory 
is approximately valid e.g. for bosonic 85 Rb atoms con- 
trolled by a Feshbach resonance [la ]. However, in gen- 
eral, bosonic quantum evaporation and three-body com- 
bination are best accommodated by generalising the GP 
equation to include non-local terms and terms of higher 
order respectively (l7j . which interfere with the simple 
form of (fT]). With these qualifications experiments on 
particle production can be performed whose outcomes 
can then be predicted (e.g. as in [l7j ) but they do not 
sit comfortably in an analogue gravity framework. 

This problem at least is largely solved for cold Fermi 
gases, for which the Pauli exclusion principle strongly 
suppresses three-body effects [l5j], permitting robust 
mean- field approximations [l8j], and it is these that we 
shall consider here. Some qualification is necessary, since 
a semiclassical representation for cold fermionic conden- 
sates is, of itself, not enough to lead to simple geodesic 
equations (fT]) for the phonon. As we have shown else- 
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where [19j |. the hydrodynamic approximation for a con- 
densate of paired cold Fermi atoms is, in general, a two- 
fluid model which does not permit a representation like 
([T|) with a single metric. A s ing le-fluid model requires a 
dominant narrow resonance |20(, such as the resonance 
in 6 Li at 543. 25G, which we henceforth assume. This 
gives us the additional benefit that, the narrower the res- 
onance, the more reliable is the mean-field approximation 
that we shall use [lH ■ 

It is also the case that, for cold Fermi gases, the linear 
behaviour of the long-wavelength condensate dispersion 
relation can also be derived from a GP equation under 
some circumstances 120 ] . Although identical to the equa- 
tion invoked in [l2T - [l4| (in terms of the speed of sound) , 
the short-wavelength non-linear behaviour, crucial for 
understanding the applicability of (JT|), now has a very 
different origin. As we shall show, the speed of sound is 
changed through the coupling of the gapless phonon field 
to density fluctuations, the collective mode with a finite 
gap [2l[. The deviation from linearity in the dispersion 
relation is due to these new degrees of freedom in the sys- 
tem, reflected in the presence of a Higgs field (22|. This 
is very different from the way that the non-linear com- 
ponent of the dispersion relation arises as a consequence 
of the quantum pressure in conventional Gross-Pitacvskii 
theory, for which the only degree of freedom is the gapless 
Goldstonc phonon. 

The plan of this paper is as follows. In the next sections 
we show how the geodesic equation ([T]) can be derived 
for low energy-momentum phonons in cold Fermi gases 
controlled by a narrow resonance. In the remainder of 
the paper we examine to what extent this equation can 
be justified for modelling phonon production as the speed 
of sound is varied. Our conclusion is that ([T]) is deceptive 
in practice because of the hidden length and time scales, 
not present in it, that are a consequence of the underlying 
Galilean invariance of the system. 



II. CONDENSATES FROM COLD FERMI 
GASES 

We begin as in our previous papers 0, [HI • At 
temperature T = in the narrow resonance limit the 
cold Fermi gas is described by the action (in units in 
which h = I) 



binding energy 2fi — v which determines the strength of 
the interactions, and which is controlled by an external 
magnetic field. Weak fermionic pairing gives a BCS the- 
ory of Cooper pairs, and strong fermionic pairing gives 
a BEC theory of diatomic molecules. The crossover is 
characterised by the divergence of the s-wave scatter- 



ing length as oc (2/x 



}18l ] as it changes sign [2f 



On driving the condensate from the deep BCS regime 
(as < 0) to the deep BEC regime (as > 0) by ramp- 
ing the external magnetic field H the speed of sound c 
decreases from 0(v F ) to essentially zero. 

Integrating out the quadratic Fermi fields gives us an 
exact non-local one-loop effective action Snl in terms of 
4>(x) = - \4>(x)\ e i6 ^ [lHH. The dynamics is encoded 
in 9. The action possesses a U(l) invariance under 9 — > 
9 + const., which is spontaneously broken by spacetime 
constant gap solutions |</>(:r)| = \(/)q\ ^= 0. 

Wc restrict ourselves to the mean-field approximation, 
the general solution to SSnl = 0, valid for a suffi- 
ciently narrow resonance [TH, [2(| . The Galilean invari- 
ants of the theory are the density fluctuation 6\<fi\ = 
= 9+(V9) 2 /4m, a ndX(5\<f>l9) = (5\<f>\) + 
V9.V(S\cj)\)/2m. 9(x) is not small. Expanding Snl in 
powers of them gives [HI, [24j a local Galilean invariant 
effective Lagrangian density L e ff for long- wavelength, 
low-frequency phenomena of the form 
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J eff 



Po G(9,e) 



-aeG(0,e) + -7A 2 (e, 



MV, (3) 



in which e oc 8\(j)\ i.e. is a density fluctuation. The scaling 
is chosen so that, on extending G(9) to G(9, e) = 9 + 
(V6*) 2 /4to + (Ve) 2 /4m, e has the same coefficients as 9 
in its spatial derivatives. For what follows, the details of 
the scaling are immaterial. In ([3]) Nq is the density of 
states at the Fermi surface and po is the total fermion 
number density. The definitions of a,j,M can be found 
in [HIHIl- The overbars denote renormalised quantities. 
Rather than repeat them here we give exemplary plots of 
the relevant quantities in Fig.l, to be discussed in more 
detail later. 



A. The hydrodynamic limit 
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The ipo- denote the Fermi fields with spin label a = (t, I). 

The diatomic field <f> describes the narrow bound-state 
(Feshbach) resonance with mass M — 2m and tunable 



The Euler-Lagrange (EL) equation for 9 is the conti- 
nuity equation of a single fluid, 



d 

—pe + V(p e v ) = 0, 



(4) 



where pe = po + 2ae — N G(9) and v# = V#/2m. 

The Euler-Lagrange equation for e is less transparent, 
of the form, 



j t [r]X{e, 9)] - £-V 2 e + M 2 e + 2aG(9) = 0. (5) 
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FIG. 1: The curves show a 2 (red), N (grey) and a 2 / M 2 
(olive) for the value g = 0.9 (defined later in terms of g) as a 
function of 1/kFO-s- 



However, if we now neglect both the spatial and tempo- 
ral variation of e, in comparison to e itself (the hydro- 
dynamic, or acoustic, approximation) the e EL equation 
becomes 

e«-2aG(fl)/M 2 . (6) 
Inserting this in ([S]) leads to the Bernoulli equation 

1 



mwg + V 



Sh H — mvg 



0, 



(7) 



where the enthalpy is Sh — 5p/ p = mc 2 8pj p. The result- 
ing equation of state is dp /dp = mc 2 across the whole 
regime. 



B. The acoustic metric 

If we now insert ([6]) into the continuity equation ((4]) it 
can be rearranged into the form 



d 
dt 



Po 
2mc 2 



(G(9)V9) 



Pove 



where 



po/2m 
N + 4a 2 /M 2 ' 



(8) 



(9) 



If we linearise (|8j with regard to 9 in a time-independent 
condensate it takes the form 



V 2 6> = 0, 



(10) 



with an effective Lorentz metric, enabling us to identify 
c as the speed of sound. 

More generally, let us suppose that the condensate has 
a background velocity Vo = V6*o/2m in terms of a back- 
ground phase 9q. HO is the fluctuating phase around the 
background, 



9 = 9 Q + 9, 



(11) 



then 



G{9) = G(9 ) +9 + v o .V0 = G(9 ) + X(6, 9), (12) 

in terms of X given earlier. 

If we now insert (fT2|) into ([5J and linearise with respect 
to 9, ([8|| then becomes, on using {7}, 



d_ 

dt 



Po 



X( 



§(v x(M)) 



PoV9 



(13) 

On rearrangement this can be written in the compact 
covariant form 



d^dj) = 0, 



in which 

r 



f« . fOl _ _ ru I . flU _ i 

' c 2 ~~ c 2 ~~ c 2 



p Q 5 lJ - 



Po i j 



We can rewrite (|T4|) as the geodesic equation 
1 



-MV^g 9^9,6) = 0, 



-9 



(14) 



(15) 



(16) 



with which we began, where in d spatial dimensions 

This is the canonical form presented in @, [l2| ■ 

C. Gross-Pitaevskii equations 

In what follows we shall restrict ourselves to static 
homogeneous condensates with constant po for which 
vo = 0. As we noted elsewhere (24|, for such condensates 
these hydrodynamic results can be derived from a Gross- 
Pitaevskii (GP) equation on ignoring quantum pressure. 
Consider the Lagrangian describing the wave-function ip 
of a particle of mass 2m, interacting non-linearly with 
itself, 

*2 2 
• n td f 

L(V0 = iff - — ■ - — (H 2 - Po) 2 (18) 
4m po 

where we have restored factors of h. The Gross-Pitaevskii 
equation following from ([T5]) is 



ihip H V 2 ifj + 2mc 2 ip ■ 

2m 



2mc^ 
Po 



0. 



(19) 
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If we set ip = yfp cxp(i9) and solve (fl~9|) at the relevant 
order in derivatives, we recover (j4]) and ((7]). We stress 
that it is through the slaving of the gapped mode (density 
fluctuations) e to the gapless mode 9 that we can describe 
the system ([3]) by a GP equation with its fewer degrees 
of freedom. 

On introducing a time-dependent phase in ip to elimi- 
nate the linear term in (j!9|) this GP equation is formally 
identical with the GP equation proposed for tunable con- 
densates in [l2l - [l4l ] , with Lagrangian density 



L'(ip) =ihip*i}j 



4m 



Wip* ■ Vip 



mc 



A) 



W<*V0 2 - (20) 



The immediate difference between the two GP equations 
lies in the relationship between c 2 and as- For the case 
of L' , corresponding to a system of bosonic atoms (l2| - [l4l | 
c 2 is linear is as as 

mc 2 Anh 2 as 



2m 



(21) 



In our case of fermionic atoms, the typical behaviour 
of c 2 as a function of as has the very different behaviour 
of Fig. 2. In Fig.l we show a 2 ,N and a 2 /M 2 for ex- 
emplary values of the condensate parameters, where we 
have renormalised the coupling strength g of ([3]) to g as 
defined later. In the deep BEC regime, where Nq is small 
and 4a 2 /M 2 large, c(t) vanishes. To a good approxima- 
tion in, 

c 2 {a s ) = (4 CS /2)[1 + tanh(d(<?) - b(g)/k F a s )}, (22) 

over the whole range from deep BCS to deep BEC be- 
haviour for a wide spread of couplings g and it is this 
parametrisation that we shall use to motivate our results. 



D. Tuning the condensate 

It will be convenient to rewrite c 2 as 




FIG. 2: The dotted line shows c 2 for the value g = 0.9 as a 
function of l/fc^as. The solid line shows the parametrisation 
(122)) for d(g) = -1.15 and b(g) = 0.202. We get as good or 
better fits for other values of g, with b(g) varying by only 25% 
over the range 0.2 < g < 1.6. 



changes the s-wave scattering length as as 2(i — v = 
g 2 No/kpas- In turn, it is determined by the value of 
the applied external field H. as fi~8j 



as = dbg I 1 - 



H — Hr 



(26) 



where a^g is the background (off-resonance) scattering 
length and the so-called 'resonance width'. Ho is the 
field required to achieve the unitary limit (\as\ — > 00). 
For the case of interest we pass from the BCS to the 
BEC regimes as H decreases through "Ho- 
For our narrow resonance 1181 



c 2 = Be 2 



BCS- 



(23) 



in m 



cbcs 



^Jpo^mN^s = vp /^ (24) 



is the velocity in the deep BCS regime (vp is the Fermi 
velocity) and 



B = 



N BCS 



N + 4a 2 /M 2 



(25) 



Nq 80 is the value of N in the deep BCS regime (see 
Fig.l). 

Applying an external magnetic field T-L changes c by 
changing the binding energy 2fX — v which, in turn, 



1 



1 



kFds kpabgWi 



-(H — "Ho), 



(27) 



If we adopt (|22j) the resulting equations can then be 
solved semi- analytically. We have seen in [2(| that, on 
comparing (|22[) to the numerical values obtained from 
(J9j) , it provides behaviour that is better than merely qual- 
itative, even though the tan/i-behaviour overestimates 
the rate at which the speed of sound vanishes. This sug- 
gests that (|2"2l will overestimate phonon production. 



III. MIMICKING FRW UNIVERSES 

Consider a homogeneous condensate in a homogeneous 
magnetic field H(t) varying in time. The phonon field 9 
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then satisfies ([T]) with the metric g(t) of (jl~7|) . in which 
c 2 (t) of is controlled by H of ([27]) . To make the situ- 
ation simpler from the viewpoint of FRW analogue grav- 
ity, we follow the authors of fl2| - [l4| in assuming that the 
system is essentially two-dimensional (a pancake conden- 
sate) . This permits us to use the given speed of sound in 
a 2D setting with a direct correspondence with the FRW 
metric. A more general analysis would throw up the same 
generic results, but with technical complications. 
As in [l2|, taking d = 2 in {UJ gives 



9nu = 
If we write 

then a 2 (t) = B{t)- 



Po 

cbcs 



-BCS 
"1^2 





a(t) 2 5ij. 



-BCS 



B(t), 



(28) 



(29) 



A. Constant quench rate - expanding universe 

We first consider one of the simplest quenches to im- 
plement experimentally in which H decreases uniformly 
in time as H/H\u = ~~ th 1 - The time dependence of 
B(t) as H changes is given from (|2"Tj) as 



B(t) = (1/2) [1 + tanh(d( 5 ) - t/ TQ )}. 



(30) 



t = is the time at which the system is at the unitary 
limit, and 



b(g)H 



(31) 



The quench parameters are related to the width of the 
resonance T by [l8[ T f» 4m/x 2 3 a 2 ff 'H 2 J /?i 2 , where p B is 
the Bohr magneton. In practice, it is more convenient 
to work with the dimensionless width 70 ~ y/T^/ep. If 
To = Ti/ef, the inverse Fermi energy (in units of K), the 
relaxation time of the shortest wavelength modes, then 



ill 

TO 



\x B H h 



■7o- 



(32) 



To be concrete, consider the narrow resonance in 6 Li at 
Ho = 543. 25G, discussed in some detail in (27} ■ [This is 
to be distinguished from the very broad Fcshbach reso- 
nance in 6 Li at 850 G] As our benchmark we take the 
achievable number density po « 3 x 10 12 cm -3 , whence 
t F fa 7 x lO~ 11 eV and 70 ~ 0.2. In terms of the dimen- 
sionless coupling g, where g 2 = (QAt 2 F / 2>kp)g 2 , 6 Li at the 
density above corresponds to g 2 J$ 1. In practice, 6(g) is 
very insensitive to g, varying by less than 25% over the 
range g = 0.1 to g — 0.9. For a condensate of density p 
it follows that g 2 cx (po/p) 1 ^ 3 and 



in 

Til 



n\po 



(33) 



where T-L is measured in units of Gauss (ms) _1 . Experi- 
mentally, it is possible to achieve quench rates as fast as 
ii ?s Q.lG/ms [27j . As we shall see below, we need such 
fast quench rates if we are to shake off Galilean invari- 
ance. To make calculations simple, we take p = po in the 
discussion that follows. 

The FRW scale factor a(t) is, from the above, 



a 2 (t) = B(t)~ 



1 



,-2d 2t/TQ 



(34) 



We have normalised a(t) to unity in the deep BCS regime. 
Since d f=a —1.15 then, insofar as the hydrodynamic ap- 
proximation is reliable, in the BEC regime (t > 0) we 
have 



o(t) 



-d e t/T Q 



(35) 



corresponding to an effective de Sitter universe with Hub- 
ble parameter H = 1/tq. That is, the simplest experi- 
mental situation of a constant quench rate looks to give 
one of the most interesting analogue models! 

The authors of fl2l[l4| also considered instantaneous 
changes in sound speed, which in our case would formally 
correspond to replacing the smooth curve of Fig. 2 by a 
step function. Since tq = O(tq) is the dc facto definition 
of an instantaneous quench, experimentally we are only 
an order of magnitude away from effective instantaneous 
change (and could improve upon this by changing the 
density). 

As an aside we note that, for cold bosonic gases, such as 
85 Rb the same parametrisation (|2l)|) of scattering length 
against external field holds (although c(t) oc as). Even 
if it were the case that the GP equation ([2U]) is valid, a 
stable condensate can only exist when as > 0, i.e. for 
Hq < H < Ho + Hu in the case of negative ab g as here. 
To try to tune H(t) to produce a de Sitter- like metric 
without causing collapse of the condensate must surely 
be impossible. 

However, from the viewpoint of analogue modelling 
we are not committed to mimicking expanding universes. 
Contracting or oscillating 'universes' are equally accept- 
able in principle, although we shall see that there are 
difficulties. 



B. Oscillating fields - cyclic universes 

We conclude this section with a brief comment on in- 
ducing an oscillation in the speed of sound through an 
oscillatin g ex ternal field, leading to the metric of a cyclic 
universe [l2j . This is an interesting case in that, with 
the final and initial states identical, particle production 
occurs only as a result of parametric excitation. 

For reasons that will become clear later, we wish to 
stay within the BEC regime (as > 0) without extending 
too far within it. There is considerable freedom in how 
we do this, but a simple choice is to take 



1 



1 



k F a s 



ip a s 



q-(1 — sinwt), 



(36) 
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corresponding from (|27|) to an oscillating field in which 
the oscillation of a$ 1 about (a^) -1 extends to the unitary 
regime. 
As a result 



a 2 (t) = 1 + a 2 



-2A sin tot 



a n e 



-2A sin cut 



(37) 



where A = b/kpas and a 2 , = exp — 2(d — b/kpas)- With 
d = —1.15 from Fig. 2 the approximation is justified. 



IV. THE BREAKDOWN OF THE SIMPLE 
ACOUSTIC MODEL 



regime, for which a(t) varies from approximate constancy 
to approximate linearity in t, in comparison to the expo- 
nential growth of (|35|) . For this and other reasons that 
will become clear later, in the first instance we restrict 
ourselves to the BEC regime (as > 0, t > 0) and for 
which Nq < 4<5 2 /M 2 (see Figs. 1). In this regime a 
is approximately constant. We find that the simple ap- 
proximations that we shall make are good enough for 
establishing more than qualitative behaviour. We stress 
that, unlike the case of cold bosonic atoms, we are not 
restricting ourselves to as > for reasons of stability, 
since condensates of cold Fermi atoms do not collapse 
for as < 0. 



We shall not attempt to perform any calculation of 
phonon production with respect to the metrics above 
since, in practice, this simple picture is never imple- 
mented. These hydrodynamic equations have arisen from 
the non-linear EL equations by ignoring derivatives of the 
density fluctuations. On the other hand, for small fluctu- 
ations, the linear approximation (while retaining deriva- 
tives of density fluctuations) to the EL equations for 9 
and e following from (|3]) is 



— — — - — V o — ae 

2 4m 

2 4m 2 



(38) 



where, for simplicity, we have assumed constant coeffi- 
cients. The residual Galilean invariance lies in the single 
time derivatives. The identification of the density fluctu- 
ations e with a gapped mode is clear from (|3"8"|) , although 
we shall take this no further. A proper discussion of this 
mode demands that we take the two-fcrmion cuts into 
account @, [H, [H| . 

On diagonalising (f3"5)) , we see that for long wavelengths 
the phonon has dispersion relation u 2 = c 2 k 2 for the 
identical speed of sound c as in ([§]) , despite the presence of 
derivatives of e that are absent in the hydrodynamic ap- 
proximation (which corresponds here to e = —2a9/M 2 ). 
Yet again a more careful analysis shows that this result 
does not require constant coefficients. 

We see that, without having to adopt the hydrody- 
namic approximation, tuning the speed of sound is im- 
plemented microscopically by the coupling of phase (ve- 
locity) fluctuations to density fluctuations. Because of 
this, the intrinsic Galilean invariance introduces time 
and length scales not visible in the 'Lorentzian' limit of 
([T]), but there nonetheless. We now consider the circum- 
stances in which they cannot be ignored, spoiling the 
emergence of an FRW geometry. This is a general prob- 
lem, that has been discussed in part (but then more gen- 
erally) in [H E3, but our model allows us to address 
these issues in a way that is wholly determined from the 
specific microscopic dynamics of cold Fermi gases. 

We first note that phonon creation would be expected 
to be less important in the BCS and the intermediate 



A. Rainbow metrics 

The underlying Galilean invariance imposes deviations 
from a linear phonon dispersion relation, most simply 
expressed as a variation in c with wavelength, once we 
get away from the long wavelength limit. Most simply, 
we follow [l3|, [HI in introducing a wavelength dependent 
metric termed the 'rainbow' metric. 

To see the dependence of the speed of sound on wave- 
length we ignore the second time derivative of e in (|38|) 
as before, but retain the spatial derivatives, replacing the 
equations on a mode by mode basis with 



No 
2 

4m 



V a — ae = 

4m 

e + -M 2 e + a9 



I) 



(39) 



for the mode with wavelength k. 

Eliminating e then gives a mode-dependent speed of 
sound 
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p /2m(l + p Q k 2 /2mM ) 



N (l + p k 2 /2mM ) + Aa 2 /M 2 
which gives the expansion around c (= Cfe = o) as 

k 2 4a 2 c 2 



4 



i 



M 2 M 2 



In the BEC regime this can be written as 



4 



Po 
2m 



4a 2 



1 



k 



2 1 



PO 

2m) Aa 2 c 2 



(40) 
(41) 

(42) 
(43) 



Away from the unitary limit in the BEC regime, a 
is approximately constant. After taking the scaling of 
e into account, we find a ~ po. Provided c is not too 
small, when (|43[) breaks down, this gives us a transitional 
momentum 



2ac 
p 0/ 1 2m 



4mc 



(44) 
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above which the short wavelength modes see the non- 
linear effects of the dispersion relation. 

To see the conse quen ces of this in more detail we repeat 
the analysis of (13, Il4| . For modes of wavenumber k the 
FRW scale a(t) is modified in a mode dependent way 
from (EH) to 



for K above. 



a 2 k (t) « [B{t) + k 2 /K 2 



B. The linear quench 



(45) 



In the BEC regime we find 



z 2d e -2t/r Q + K r -2d c 2t/To 



K 2 



where 

K Q = 4macBcs/po ~ 4mc BCS = (4/V3)k F , 



(46) 



(47) 



with kp the Fermi momentum. From our earlier com- 
ments, we assume that t > 0. Some caution is required 
in that the approximation (|46| , derived from (|41j) , breaks 
down when the second term is too large. 

Scale factors which may begin as expanding can stall 
and rapidly contract, as is seen from the 'rainbow' Hubble 
parameters in the BEC regime, 



H k {t) = 



Qfc(t) M ^1 

Ofc(t) Tq 



g2dg — 24/tq _j_ k ^ 



-16. It/TQ 



z 2d e -2t/r Q _|_ J^ e -2d e 2t/TQ 



(48) 

Taken literally, Hkit) decreases monotonically, flipping 
rapidly from positive to negative values. The transi- 
tion between an expanding and collapsing 'universe' for 
wavenumber k happens at time tk, 



tk 



d + In 



k ' 



in terms of which 
H k (t) 



-1 



tanh(2 (t - tfe)/rg). 



(49) 



(50) 



From our comments above this approximation breaks 
down for t>t k . 

There is an infrared bound on k. For the exemplary 
condensates with N « 10 5 atoms, their width is £o ~ 
10 2 /fc F 0. This gives 



0< — <d 



In 10, 



or perhaps a little larger. With d w — 1 this gives 
0<t fc <r Q , 



(51) 



(52) 



at best (without having to restrict ourselves to the BEC 
regime a priori). Taking the example of Fig. 2 with b ~ 
0.2 this translates into a transition between an expanding 
and a collapsing 'universe' for a value of t for which 



< l/a s k F < 5, 



(53) 



as we go from the shortest to the longest wavelengths. 
We stress that tk (or the corresponding 1/asfcp) marks 
the boundary between the applicability of Eq.([T]) with 
its 'Lorentzian' structure and the restoration of Galilean 
invariance. 

As a guide to phonon production during the sweep 
from BCS to BEC we also need the modified dispersion 
relation 



UJ k (t) 



e 2d e -2t/r Q 



K 2 



e -2d e 2t/r Q 



1/2 



V2u e d e-* fc / T «{cosh(2(t-^.)/r Q )} 1 /^54) 



where we have adopted the notation of [l3|, [l4| , in which 
luq = \k\cBCS- The relevant quantity is the ratio 



n k {t) 



H k (t) 



k {cosh(2(t-t k )/T Q )} 3 / 2 



2TQCBCsk VK sinh(2(t-t fe )/r Q ) 

2 tq( k \ 3/2 {cosh(2ft-^.)/r Q )} 3 / 2 
V3 T V k F J sinh(2(t - ifc)/r Q ) 



,(55) 



using our definitions of tk and Kq given earlier. 
A quantum mode with wavenumber k only experiences 
significant amplification (and hence phonon production) 
when lZk{t) « 1. As before, this approximation breaks 
down when t > tk- 

In the vicinity of tk, where Hk is small, IZk is cor- 
responding large, as it is for t much greater than tk- 
In between it achieves a minimum at t = t*, where 
cosh{2(i* - t k )/T Q ] = V3. That is, 



t k 



-tq/2, 



at which 



n k (t* 



r V k F 



3/2 



(56) 



(57) 



[The other minimum at t* — t k ~ tq /2 is unreliable.] In 
order to have any phonon production we must have as 
fast a quench as possible, with a current lower bound of 
tq/tq w 10 and a lower bound of k/k F of 10~ 2 , say, for 
our typical condensate. Then, for the lowest momentum 
phonons, 



K k (t*) « 0.04, 



(58) 



this minimum increasing as momentum increases. Thus, 
from (|55p . there is a window in which lZ k (t) is sufficiently 
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small to expect phonon production. However, this win- 
dow is not large and we note that, if these approxima- 
tions were reliable for t > tk then, for such low momen- 
tum phonons we would expect comparable production in 
both the expanding and contracting phases. There may 
or may not be a truly contracting phase but, whatever 
the details, since the contracting phase is purely a con- 
sequence of Galilean invariancc, it is clear that we are 
not deriving the particle production that we might have 
anticipated from Eq.(fT]). For higher momentum phonons 
and somewhat slower quenches lZk{t) 3> 1 throughout 
and there is no phonon production. 

As for reversing the direction of the quench, we have 
seen that the non-linear effects are greater as c becomes 
smaller, making a quench beginning in the deep BEC 
regime problematical. 



Periodic field 



For the periodic quench discussed earlier we have 



4(t) 



-2 „2ytsinui 



n 2 -2Asinu>t 
(In £Z 



2 "0 



(59) 



from which the rainbow Hubble parameters follow as 



H k (t) 



-Auj cos u>t 



r —2 2Asinujt 



fc n 2 „-2Asinujt 
K? a e 



n ~ 2 ^Asinujt j_ fc 2 ^,2 ^ — 2 A sin ujt 

a e t -^t" e 



(60) 

not at all simple, despite its periodicity. We shall not pur- 
sue this further here, beyond noting that, from (|6T)|) . the 
transitional momentum marking the boundary between 
Lorentzian and Galilean behaviour is (with A ~ 0.1) 



these linearised EL equations display the oscillatory be- 
haviour 



s p B = sp B cos sit, 



(62) 



describing the repeated dissociation of molecules into 
atom pairs and their reconversion into molecules. The 
frequency fl of density fluctuations is determined by the 
energy scale at the beginning of the two-fermion cuts 
in the energy plane (E m i n of [28|). It can be shown 
[28|,[29| that fi increases monotonically from the exponen- 
tially damped 0(/xexp(— Tr/2kp\as\) in the BCS regime 
to O ps 2/i in the deep BEC regime ft w 2/i. For 
the single-fluid model to be valid, the density fluctua- 
tions must be able to be averaged to zero on the natural 
timescale r = fi/Mc 2 of (TTg|) i.e. rfl ^> 1. This is read- 
ily achieved in the BEC regime, but is difficult, if not 
possible, to achieve in the BCS regime [2(|. For our ex- 
emplary condensate we find tQ « 2tt when l/kpo-s ~ 0, 
falling fast as 1/kpas goes negative e.g. rfl « 1 when 
1/kpas ~ — 1 (achieved at time t « — tq/5 for the fastest 
quenches). 

We note that a further reason for our restriction to t > 
is that, in the discussion above, we have taken a ~ po 
constant and large. As can be seen from Fig.f, a falls 
away near the unitarity point. Even had it been sensible 
to take t < 0, reducing a reduces Kq, effectively making 
low momentum phonons behave as UV phonons, reducing 
or eliminating their production in the BCS regime, as we 
had already anticipated. 

In this regard at least, there is no further problem with 
oscillatory universes in the BEC regime, for which we 
would require ft ^ w. This is automatically the case 
for the system considered here, since simple calculation 
shows that tqCI >• 1 for the BEC regime, and we must 



have 



> t . However, this needs further study. 



n 2d 



K 



10' 



(61) 



or somewhat less. With 10~ 2 < k/Ko < 1 for our exam- 
ple this means again that the bulk of the activity is con- 
trolled by the Galilean group, a conclusion aided by our 
earlier observation that the parametrisation (|22[) overes- 
timates particle production. 



D. Time scales 

For the linear quench we have already noted that, as we 
move to the BCS regime and away from exponential be- 
haviour, iJfc diminishes, making phonon production less 
likely. However, the reason why we have to restrict our- 
selves to the BEC regime (i.e. t > 0) concerns the time 
scales hidden by the Galilean nature of the theory. The 
coupled equations (|38|) describe a two-component system 
of molecules and atom pairs. The density of molecules is 
p B = 2|0| 2 , whose fluctuation is Sp B oc e. For our spa- 
tially homogeneous condensate, ignoring damping [23[, 



V. CONCLUSIONS 

We find the idea of the emergence of relativistic 'space- 
time' from condensed matter systems, as manifested in 
the idealised Eq.([T|), beguiling. However, in this paper 
we have shown that the underlying Galilean invariance 
of realistic condensed matter systems can impose a brutal 
reality check. 

We have explored this in the context of bosonic con- 
densates of cold Fermi atoms tunable through a nar- 
row Feshbach resonance (a broad resonance requires a 
doubled metric). To be specific we have considered 6 Li 
atoms in concentrations and volumes that arc experimen- 
tally accessible and with quenches that are equally well 
achievable with current practice. To have an analytic 
'emergence' requires that we have a robust analytic semi- 
classical approximation and cold Fermi gases permit this 
in a way that, in general, bosonic atoms do not seem to 
possess because of many-body effects. In addition, the 
condensate degrees of freedom for Fermi gases, with a 
gapped 'Higgs' mode, arc very different from those of the 
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elementary Bose condensates invoked in [12l - H4| . which 
have provided the counterpoint to this paper. It is only 
insofar as the gapped mode is slaved to the phonon mode 
that the Fermi system can be represented by the same 
Gross-Pitaevskii field as the bose system, with its limited 
degrees of freedom. This hides the fact that the speed 
of sound in condensates of Fermi atoms changes through 
the coupling of the gapless phonon field to density fluctu- 
ations, the collective mode with a finite gap. As a result, 
cold Fermi and bose gases have very different mechanisms 
for acquiring the non-linearities in the phonon disper- 
sion relation, which in the latter case is due to quantum 
pressure but, in both cases signals a restoration of the 
underlying Galilean invariance. 

From the viewpoint of analogue gravity there is no 
compulsion to mimic realistic universes; expanding, con- 
tracting and oscillatory metrics are equally acceptable. 
What matters is that we can cast them in the form of (jTJ) , 
make predictions and perform experiments to test those 
predictions. Unfortunately, once we have taken the addi- 
tional degrees of freedom of a Fermi system into account 
we have shown that attempts to mimic particle produc- 
tion in FRW universes by tuning the speed of sound are 
complicated because of the way in which the underly- 
ing Galilean invariance is embedded in the equations of 
motion. 

Our major analysis has been for the simplest experi- 
mental field quench, one of constant rate, which induces 
an acoustic de Sitter metric in the analogue FRW sys- 
tem, one of the better studied systems 12-[l4|. For low 
momentum phonons and superfast linear quenches, the 
best that can happen is that the true effective metric can 
show a shortlived de Sitter-like expansion before stalling, 
and perhaps contracting. In fact, the problem is worse 
since, for sweeps from the deep BCS to deep BEC regions 
attempting to mimic expanding FRW universes, the hy- 
drodynamic approximation breaks down before reaching 
the BEC regime, again due to the underlying Galilean in- 
variance. As a result, at most a fraction of the expansion 
may correspond to the geodesic picture of (fTJ . There are 
similar problems with oscillating fields and their analogue 
cyclic universes. There is the caveat that all the above 
is, in detail, for one 6 Li configuration and it is possible 
that very different parameter choices (including different 
Fermi gases) will give better results, despite the generic 
nature of the problem. 

The alternative of looking to tunable cold bosonic gases 
has its own problems. This is because the inclusion of 



many-body effects in a semi-classical description is still 
not fully resolved. While we should not preclude the 
possibility of a robust mean-field description in some con- 
texts, the standard way to include these effects takes us 
out of the simple geodesic description of ffl. If these 
problems can be ignored the results of [12h14| ] do give a 
narrow window for success that has evaded us, but more 
work on the robustness of the approximations needs to 
be performed. 

This should be contrasted to sweeps from the deep 
BEC towards the BCS regime, for which we have argued 
[2(i| that, rather than particle creation, spontaneous vor- 
ticity can appear as a consequence of causal horizons. 
This relies on causality alone, without the need to con- 
struct an effective metric. 

Finally, we have to say that, despite our attempts to 
avoid idealisation of condensates and to introduce the 
complexities which, through Galilean invariance, wreck 
the analogue gravity programme, we are still dealing with 
oversimplified systems. Although the early universe may 
have been homogeneous, condensates are not, requiring 
traps to contain them. We have avoided discussing traps 
on the grounds that, if we cannot get a homogeneous 
system to do what we would like, there is little need to 
worry about further complications. 

We stress that these concerns do not apply to other 
work by several of the same authors of |12Hl4j . in 
which they propose physically realisable highly tuned 
two-component condensates pOl . l3~l| as a way to create 
massive degrees of freedom, rather than through the 
Higgs mechanism, as here. This may be a productive 
way to proceed. 
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